Remarks on possible local parity 
violation in heavy ion collisions 

Adam Bzdak a ' 6 , Volker Koch a , and Jinfeng Liao° 

Lawrence Berkeley National Laboratory, 1 Cyclotron Road 

MS70R0319, Berkeley, CA 94720, USA* 
b Institute of Nuclear Physics, Polish Academy of Sciences 
Radzikowskiego 152, 31-342 Krakow, Poland 

December 27, 2009 



Abstract 

In this note we discuss some observations concerning the possible local 
parity violation in heavy ion collisions recently announced by the STAR Col- 
laboration. Our results can be summarized as follows (i) the measured corre- 
lations for same charge pairs are mainly in-plain and not out of plane, (ii) if 
there is a parity violating component it is large and, surprisingly, of the same 
magnitude as the background, and (iii) the observed dependence of the signal 
on the transverse momentum (pt) is consistent with a soft boost in pt and 
thus in line with expectations from the proposed chiral magnetic effect. 
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1 Introduction 



Recently the STAR collaboration announced pQ the results on possible local parity 
violation in heavy ion collisions. In Ref. [H [3] it was argued that in the hot dense 
matter created in heavy ion collisions local, instanton or sphaleron, transitions to 
QCD vacua with different topological charge may result in metastable domains, 
where parity is locally violated. 

In this paper we will solely concentrate on an analysis of the experimental re- 
sults. We will neither attempt to provide alternative explanations for the observed 
correlations, such as e.g. given in Ref. [I] nor will we discuss the likelihood that the 
proposed effect may occur in a heavy ion collision. For detailed discussion of the 
underlying mechanism and the latest theoretical review of this problem we refer the 
reader to Ref. [3]. 

The phenomenon due to local parity violation, which is of relevance for the 
discussion here, is the so called chiral magnetic effect [21 [3]. It leads to the separation 
of negatively and positively charged particles along the system's angular momentum 
(or equivalently the direction of the magnetic field) into two hemispheres separated 
by the reaction plane. As a result, the system exhibits a electric current along the 
direction of the angular momentum, and thus breaks parity locally in a given event. 
However, since instanton (sphaleron) and anti-instanton (anti-sphaleron) transitions 
occur equally likely, the chiral magnetic current is either aligned or anti-aligned 
with the angular momentum. As a result, the expectation value of any parity odd 
observable, such as \JomI^ vanishes. Here jcM is the chiral-magnetic current and 

I is the angular momentum. Consequently, a direct measurement of parity violation 
even in a small subsystem is impossible. However, one may attempt to identify 
the existence of these parity violation domains by studying the fluctuations or the 
variance of a parity-odd observable. Since the variance of a parity-odd observable 
is parity even, in principle other, genuinely parity-even, effects may contribute, and 
one needs to separate those carefully before being able to draw any conclusions 
about the existence of local, parity violating domains. 

In Ref. [5] Voloshin proposed a method to measure the variance of a parity odd 
observables. He suggested to measure the following correlator (cos(0 a + (pp — 2^/ RP )), 
where ^rp, <ft a and 0^ denote the azimuthal angles of the reaction plane and pro- 
duced charged particles respectively, see Fig. [fl 

As we will discuss in more detail in Section [2] this rather involved correlation 
function has the advantage that correlations which are independent of the reaction 
plane do not contribute. As a result, a large fraction of the expected background 
should cancel. Recently the STAR collaboration has reported the measurement of 
the above correlation function [TJ , both integrated over the entire acceptance as well 
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Figure 1: The transverse plain in a collision of two heavy ions, ^rp, 4> a and 0^ 
denote the azimuthal angles of the reaction plane and produced charged particles, 
respectively. 

as differential in transverse momentum and pseudo-rapidity. 

This paper is organized as follows. In the following section we will analyze 
the integrated STAR result and will suggest additional measurements necessary to 
further clarify the situation. In the subsequent section, we will concentrate on the p t 
differential results and explore to which extend they are consistent with the expected 
soft phenomena due to the chiral magnetic effect. 

2 The integrated signal 

In Ref. p] the details of the STAR measurement are given. Among other things 
STAR shows the results for (cos(0 Q — <pp)) and for (cos((p a + (ftp — 20 c )), where <fr a ,p,c 
are the azimuthal angles of the produced charged particles. The paper gives reason- 
able arguments that 

(COS(0 Q + (j)f3 - 20 c )) = (COS(0 Q + 0/3 - 2^ RP )) V 2>c , (1) 

where ^rp is the angle of the reaction plane, and t>2, c characterizes the elliptic 
anisotropy for the particle with angle C . 

For the rest of the discussion we will assume that the relation ([1]) is correct. As 
a consequence, we will work in a frame where the reaction plane is defined by the 
x — z coordinates and where the y direction is perpendicular to the reaction plane. 
In other word we work in a frame where ^ rp = 0, see Fig. [U Furthermore, since 
(cos(0 Q — (f)fs)) is independent of the direction of the reaction plan^j], it will be the 
same also in the frame where the reaction plane is specified e.g., ^rp = 0. Thus 

indeed (cos(</> Q - cpp)) = (cos([0 Q - t> B ,p] - [4>p - * hp])) 
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within our frame we have to consider the following two-particle correlations: 

(cos(0 Q - = (cos(0 Q ) cosO/3)) + (sin(0 a ) sin(^)) , 

(cos(<p a + (j>p)) = (cos((p a ) costyp)) - (sin(0 a ) sin^)) . (2) 

STAR has measured both these correlation functions for same sign, (+, +) , (— , — ). 
and opposite sign, (+,—), pairs of charged particles. Qualitatively the data for 
Au + Au collisions can be characterized as follows. 

• For same sign pairs: 

(cos(0 Q + <pp)) same ^ (cos(0 Q - <pp)) same < 0. (3) 
Using Eq. (j2J) this implies 

(sin(i)sin(^)) s « - 0, 

(cos(0 Q ) cos(^)) somc < 0. (4) 

• For opposite sign pairs we find that 

(cos(0 a + fo)) opposite ^ 

(cos(4> a -^)) oppostte > 0. (5) 

Again, using Eq. (j2J), this means 

(sin(0 a ) sin(0 /3 )) opposjte ~ (cos(0 Q ) cos(0 /3 )) oppos . te > 0. (6) 

The actual data decomposed into the above components are shown in Fig. [2j 

The fact that for same charge pairs the sinus-term in Eq. (j3J (see Fig. |2J is es- 
sentially zero whereas the cosine term is finite, tells us that the observed correlations 
are actually in plane rather than out of plane. This is contrary to the expectation 
from the chiral magnetic effect, which results in same charge correlation out of plane. 
In addition, since the cosine term is negative, the in-plane correlations are stronger 
for back-to-back pairs than for small angle pairs. Second, we see that for opposite 
charge pairs the in- and out-of-plane correlations are virtually identical. This is 
hard to comprehend, given that there is a sizable elliptic flow in these collisions. At 
present, there is no simple explanation for neither of these observations. However, 
they may be explained by a cluster model, which requires several, not unreasonable, 
assumptions [1]. 

One may ask if there is room for a parity violating component if for the same 
sign (sin ((/> Q ) sin (</>£)) ame ~ 0, i.e. the signal is in-plane rather than out of plane. 
Following the argument of Ref. [HE], we can always write 

(sin(0 a ) sin(^)) sQme = B out + P, (7) 
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Figure 2: Correlations in-plane (cos(0 a ) cos^p)) and out of plane (sin(0 a ) sin^^)) 
for same and opposite charge pairs in Au + Au collisions. As can be seen the 
correlations for same charge pairs are mainly in-plane. 

where P is the part of the correlation which is caused be the parity violation (at this 
stage we do not claim that P ^ 0) and B out represents all other contributions by 
correlations projected on the direction perpendicular to the reaction plane. Denoting 
the correlations in-plane (cos(0 Q ) cos((pp)) same by B in we obtain: 

(COS(0 Q + (f)p)) same = [B in ~ B out ] - P, 

(cos(0 a - (j)p)) same = [B in + B out ] + P. (8) 

The advantage of (cos(0 Q + 4>p)) is obvious. The background is B in — B out , 
meaning that all correlations that do not depend on the reaction plane orientation 
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cancel. The STAR collaboration studied many known sources of reaction plane 
dependent correlations and all effects produce B in — B out which is much smaller 
than the observed signal. We note, however, that at present the background is 
not understood since none of the present models is able to explain the value of 

(COS(0 Q - (f>p)). 

Following the above argument, however, immediately implies that (using Eq. (j3J) 
and Eq. ©) 

P ~ -B out ~ -B in , (9) 

i.e., the parity violating effect has to be precisely of the same magnitude as all other, 
standard correlations. This relation is quite an unexpected coincidence. It means 
that the parity signal is quite strong, and consequently should also be visible in 
(cos(0 Q — 4>i3)) S ame ^ the background is well understood. 

In our view, it is mandatory to explore if the relation, Eq. (Q, is just a co- 
incidence or an indication of potential problems with the present interpretation of 
the data. To answer this question it is essential to analyze the correlation function 
(cos(4> a — (pp)) same differentially in transverse momentum and pseudo-rapidity as it 
has been already done for (cos(0 Q + 0/3)) same - Should relation (Q persist also for 
the differential correlations, one would have to conclude that the proposed parity 
violating effect is not seen in the data. 



3 Transverse momentum dependence 

The STAR collaboration has also presented [TJ the measurement of (cos(0 a + <ppj) 
in mid-central Au + Au collisions as a function of p + = {j) t ,a + Pt,p)/2 an d p~ = 
\pt,a ~ Pt,p\, where p t , a and p tj p are the absolute values of the particles momenta. 
Qualitatively the data can be characterized as follows^ 

• For same sign pairs in the range < p+,p~ < 2.2 GeV: 

(cos(0 Q + same oc p + 

(cos(0 a + <f>p)) p _ t same ^ const. (10) 

• For opposite sign pairs the signal vs p + and p_ is consistent with zero. 

One would expect [U [2] that the parity violating signal should be a soft, low p t 
phenomenon. Thus the observed increase of the signal for same sign pairs with p + 
seems to be inconsistent with the chiral magnetic effect. As we will show such a 

2 In the present Section we are only interested in the pt dependence of the signal, not in the 
overall normalization. 
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conclusion is not necessarily correct and the true signal may indeed be consistent 
with the expected low p t dynamics. 

Indeed, by definition 

(cos(0 Q + ^)) = ^p, (11) 

Wall 

where N corr is the number of correlated pairs [via cos(0 a + (pp)} and N a u is the 
number of all pairs. The latter can be easily approximated by \p + = {p t ,a +Pt„a)/2 
and p_ = \p t>a -Pt,p\)- 

Nau (p+) oc J d 2 p t . a d 2 p t ^ exp (~^r) exp (~^r) $ ( 2 P+ ~ [Pt,a + PtA) 

oc p3 e -2 P+ /T (12) 

and 

N a ii (p_) oc J d 2 p tta d 2 p tif3 exp (~^r) exp (~^r) (P- ~ |Pt,a - Pt,/9|) 

oc TV (p - /T) (p_+T), (13) 

where in the following calculations we take T = 0.22 GeV@ 

The calculated distributions of all pairs vs (p tt0l + p t ^)/2 and \p t<a — Pt,p\ are 
presented in Fig. [3j It is worth noticing that both functions are concentrated in the 
small pt region, reflecting typical thermal distributions for p_ and p + . Due to the 
soft nature of chiral magnetic effect, one expects that the distributions in p_ and 
p + for the correlated particles should not differ much from the underlying thermal 
distributions. This is indeed the case as we will demonstrate next. 

In order to estimate the distribution of correlated same sign pairs it is sufficient 
to multiply Eq. ( flOl) by the expressions ( 1T21) and ( |T3l) . respectively. Consequently 
we obtain 

N^r (P-) OC N a u (p_) , 

N CO rr{p+) OC p + N a ii(p + ). (14) 

As can be seen the dependence of the number of correlated same pairs vs 
\pt,a — Pt,p\ is identical to the dependence of all pairs presented in Fig. |3j Clearly 
the signal is concentrated in the low p t region and indeed is unchanged from a ther- 
mal distribution. In Fig. H] the dependence of the number of same sign pairs vs 
(Pt,a +Pt,f3)/2 is compared with the dependence of all pairs (previously shown in 
Fig. [3]). We find that the momenta of correlated particles are slightly shifted to the 
higher p t and the shape is roughly similar. The momentum shift required by the 



3 It corresponds to the average transverse momentum of the pions (pt) — 0.45 GeV. 
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Figure 3: The distributions of all charge pairs vs. ipt, a +Pt,p)/2 an d \Pt,a ~ Pt,p\, 
respectively. 

data is 5p + ~ 150 MeV which could conceivably be due to the large magnetic field, 
although it is somewhat on the high end of what one would naively expect from 
electromagnetic phenomena. 



4 Conclusions 

In this note we have discussed several aspects of the recent measurement of possible 
local parity violation in Au + Au collisions by the STAR Collaboration. We made 
the following three observations: 

(i) For particles with the same charge STAR sees large negative correlations in- 
plain (cos(0 a ) cos(c^g)) and very small correlations out of plain (sin(0 Q ) sin(^> ( g)) 
For opposite sign correlations in-plane and out-plain are both positive and of the 
same magnitude. 
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Figure 4: Distribution of all pairs (solid line) compared with the distribution of 
same sign correlated pairs (dashed line). Both functions are concentrated in the low 
p t region. 

(ii) If there is indeed a parity violating component in the STAR data it has to be 
of the same magnitude as all other, "trivial' correlations projected on the direction 
perpendicular to the reaction plane. This may be a pure coincidence or an indica- 
tion that the present interpretation of the data as a signal for local parity violation 
needs to be revised. To investigate this problem in more detail we need differen- 
tial distribution (vs pseudo-rapidity or transverse momenta) of (cos(0q, + c^g)) and 
(cos(0 a — (pp)) at the same time. 

(iii) We have also argued that the distribution of the number of correlated pairs 
is concentrated in the low pt region i.e. pt < 1 GeV. It is not inconsistent with the 
predictions of the chiral magnetic effect. 

At present, the data from the STAR Collaboration does not allow for a definitive 
conclusion about the presence of local parity violation. The measurement of the 
correlation function (cos(0 a — <f>p)) differential in transverse momentum and pseudo- 
rapidity is absolutely essential to further distinguish between trivial correlations and 
those due to the chiral magnetic effect. 
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